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Abstract 

Although several theories relate the steep slowdown of glass formers to increasing spatial corre- 
lations of some sort, standard static correlation functions show no evidence for this. We present 
results that reveal for the first time a qualitative thermodynamic difference between the high 
temperature and deeply supercooled equilibrium glass-forming liquid: the influence of boundary 
conditions propagates into the bulk over larger and larger lengthscales upon cooling, and, as this 
static correlation length grows, the influence decays nonexponentially. Increasingly long-range sus- 
ceptibility to boundary conditions is expected within the random first-order theory (RFOT) of the 
glass transition, but a quantitative account of our numerical results requires a generalization of 
RFOT where the surface tension between states fluctuates. 



Supercooled liquids show a dramatic slowdown of their dynamics upon cooling (14 decades 
increase of viscosity pQ within a narrow temperature range) without any obvious structural 
or thermodynamic change The dynamic correlation function of a supercooled liquid fur- 
thermore becomes progressively more nonexponential as the temperature is reduced. Several 
theories relate these phenomena to increasing spatial correlations of some sort [31 HJ El E] . 
However, static correlations have so far been unable to distinguish qualitatively the high 
temperature and the deeply supercooled liquids (though indirect evidences have been in- 
ferred from the specific heat [7] and the linear dielectric susceptibility [H]). Inspired by 
critical phenomena, it is natural to expect that the slowing down of the dynamics is related 
to the vicinity of a thermodynamic phase transition, where some kind of long-range order 
would set in [9j. This is the spirit of different recent theories [H El [TOl ttH [111 [I3], but 
appears at odds with others O [H], at least at first sight. In particular, the crucial physical 
mechanism at the root of random first-order theory (RFOT) |1| is the emergence of long 
range amorphous order, whose precise definition and quantitative characterisation is how- 
ever far from obvious. Dynamic heterogeneities [15] do show a growing dynamic correlation 
length accompanying the glass transition, both experimentally fLQ\ and numerically [171 and 
refs. therein]. This is certainly a first important step, but not sufficient to prune down — 
even at a qualitative level — different theories of the glass transition. In particular, it is 
not clear whether this phenomenon is due to an underlying static or purely dynamic phase 
transition. 

The approach followed here is based on the very definition of a thermodynamic phase 
transition, where the effect of boundary conditions becomes long-ranged. The problem is 
that for glasses there are no natural boundary conditions, since these should be as 'random' 
as the bulk amorphous states that they favor. A possible solution has been proposed in [10] , 
and further discussed in [9], in the context of RFOT, but the scope and some conclusions 
of the gedanken experiment proposed in [lUj are more general [TSl [TU] . Starting from a 
given equilibrium configuration, one freezes the motion of all particles outside a cavity of 
radius R and studies the thermodynamics of the mobile particles with boundary conditions 
imposed by the frozen ones. Defining a suitable overlap (inside the cavity) between the 
original equilibrated configuration and the configurations equilibrated with the boundary 
pinning field, the existence of "order" on a scale ^ implies a large overlap (or 'point-to-set' 
correlation) for i? <^ ^ and a small one for i? ^ ^. 
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Such correlation length arises in RFOT [TO], if one assumes (or finds by approximate 
calculations fIU{ 1^ ) that there are many amorphous metastable states, Af ~ exp[i?'^E] 
inside the cavity (S(T) is the configurational entropy, or complexity). The particles in the 
cavity can be either in the original state a in which the outside particles have been frozen or 
in any of the other states. The probability to fiip to a different state 7 is determined by the 
balance between the free-energy loss due to mismatch between a and 7, T{T)R^, and the 
gain from the entropic term TR'^Tj{T) (T is a generalized surface tension and 9 < d — 1). A 
crossover between a boundary dominated regime (large overlap) and an entropy dominated 
regime (small overlap) is expected for R = ^ oc (T/TS)^. Note that this length diverges 
at the Kauzmann temperature Tk where S(T^) 0. 

This gedanken experiment was realized numerically by three of us in ref. [221 This first 
study indeed suggested a growing static length, but this length was very small, and no sharp 
transition between high and low overlap was observed, in contrast with expectations based 
on RFOT. Is RFOT possibly in contradiction with numerical experiments? What is the 
physical meaning of the increasingly long range susceptibility to boundary conditions if it 
is characterized by such a small length? In order to answer these crucial questions we have 
performed new simulations that measure the local overlap at the center of the cavity. In 
fact, a major difficulty in interpreting the results of ref. [22]is that the overlap was computed 
as an average over the whole cavity. As a result, the overlap drop due to the (putative) 
change of state of the cavity is mingled with the decay of the overlap expected from a trivial 
weakening of the surface pinning field for larger spheres. This effect is indeed present even 
in the single state case [19j at all temperatures. As we shall show, analyzing the behavior of 
the local overlap at the center of the cavity yields sharper results which allow us to answer, 
at least partially, the above questions. 

I. BEHAVIOUR OF THE OVERLAP AT HIGH AND LOW TEMPERATURE 

We study a soft-sphere model [23] that we can equilibrate below the Mode Coupling 
transition temperature Tmc =0.226 [21] and for large systems (see methods). After equili- 
bration, several independent reference configurations are chosen as starting points for runs 
with all but M particles frozen. These mobile particles are confined inside a sphere of ra- 
dius R such that the inside density equals the bulk density. After the confined runs reach 
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equilibrium, the local overlap at the center qdR) is measured. To define qdR), we partition 
the simulation box in many small cubic boxes of side i, such that the probability of finding 
more than one particle in a single box is negligible. Let equal the number of particles in 
box i, then 

<lc{R) = ^ $^(n.(to) n,{to + oo)) (1) 

where the sum runs over all boxes within a small volume v at the center of the sphere, Ni is 
the number of boxes, and (. . .) means thermal average. To minimize statistical uncertainty 
without losing the local nature of qdR) we choose Ni = v/i^ = 125. Normalization is such 
that the overlap of two identical configurations is 1 on average, while for totally uncorrelated 
configurations qc = qo = = 0.062876. 

We show qc{R) — qo for several temperatures in Fig. [l] The decay becomes slower at 
lower temperature: the effect of boundary conditions propagates on larger length-scales. 
This clearly shows the growth of static order that expands into the bulk in the deeply 
supercooled phase. The most striking result, however, is that at low temperatures the 
decay is no longer the simple exponential that prevails in a standard liquid state [12]. The 
relaxation can be fit by a "compressed exponential" , 

qdR)-qo = ^exp[-{R/0'^], C > 1, • (2) 

where ( increases at low temperatures (see Fig. [2^) above its high temperature liquid value 
C = 1 (see Fig. [2|d). Larger values of C means a sharper crossover between large and small 
overlaps. The best fit parameters are given in table [Tj Interestingly, the value of the length 
scale ^ found here is significantly larger than those in [22] • This is partly due to the fact 
that C > 1 at low temperatures (see appendix for further discussion). The overlap qdR) 
therefore appears as a thermodynamic quantity able to single out the deeply supercooled 
liquid in terms of (a) a large correlation length and (b) an anomalous nonexponential 
relaxation, characterized by the exponent C > 1- 

The growth of C, strongly suggests that the liquid state should be thought of as a mosaic 
of regions (transiently) ordered on a lengthscale ^. Furthermore, following suggestions based 
on RFOT[ll [ini EI], it is natural to conjecture that we are probing the growth of positional 
amorphous order. Interestingly, in this scenario, all finite point static correlations remain 
featureless, whereas the point-to-set correlations [9l ES], captured by qdR)i grow. The 
simplest version of RFOT predicts at low temperature a sharp cross-over of qdR) ior R = ^ 
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FIG. 1: Overlap at the centre of the mobile cavity vs. radius R of the cavity, for temperatures 
T =0.482 (diamonds), 0.350 (triangles), 0.246 (squares) and 0.203 (circles). Lines are fits to eq. Q. 
Inset: Comparison of qc{R) — qo at T = 0.203 (filled circles) with the overlap Q{R) — qo integrated 
over the whole sphere (open circles, data ref. [22]). The local observable qdR) shows a much sharper 
behaviour. 

[26] . Our numerical results show that this is not what happens: the cross-over is much 
smoother than expected (see also the recent work [27j on a 1-d Kac model). 

II. RFOT WITH FLUCTUATING SURFACE TENSION 

To account for our numerical results, we now propose a natural generalization of RFOT, 
based on the idea that the effective interface tension, T, is in fact state-dependent. Repeating 
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FIG. 2: a) Local overlap at the centre of the mobile cavity vs. RiorT = 0.203 (circles) with fits to 
a compressed exponential (full line) and simple exponential (dotted line) . The best exponential fit 
is is rather poor. Inset: same data in semilog axes, where a simple exponential looks like a straight 
line, b) As in panel a) but for T = 0.350. The exponential fit (line) is good at this temperature, 
as can be seen in the inset, which displays the same data in semilog axes. 
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0.482 1* 0.617(40) 5.3(1.3) 
0.350 1* 0.845(28) 3.22(32) 
0.246 2.79(52) 3.04(24) 0.552(61) 
0.203 4.00(60) 3.82(12) 0.46(11) 

TABLE I: Best fit parameters as function of temperature. Errors are jacknife estimates. (*) For 
the highest temperatures, the results quoted are for a simple exponential fit, which gives a very 
good description of the data with only two parameters. If nonetheless those data are fit with the 
compressed exponential form, an exponent C compatible with 1 is obtained, albeit with very large 
(60-80%) error. 

the argument of ref . [lOl with this extra assumption, the probability that the cavity is found 
in a state 7 different from the outside pinning (frozen) state a, is: 

where /3 = 1/T, /* and S* = S(/*) are respectively the free energy and configurational 
entropy of the equilibrium states, and Pq,(T|/) is the distribution of surface tension for state 
a with other states of free energy / (see Methods). Since our simulations average over many 
realizations of the external state a, we in fact compute piut(-R) averaged over all possible 
pinning states: 

Pont[K) - 2_^Wa Poutl^J " 2^ 1 + / dT 6^"^* - Po,{T\ f*) ' ^ ' 

where Wa is the thermodynamic weight of each state. The leading contribution of the integral 
over T is shown in Methods to be e"^^"'^", where i/a is the lower edge of Pa(T\f) and where 
the exponent 6 is possibly renormalized (this might explain why the effective value of 6 can 
be less than the naively expected value d—1). Introducing P{Y) = J2a f^a^iX ^Va) we find 

d^nY).^^^^^^,^^. (5) 

The simple mosaic result of ref. [iniis recovered setting pa{T\f) = 5{T — To) and P{Y) = 
S{Y — Yq). Expression ([s]) can be simplified further by noticing that the approximation 
(1 + e^'^*-'3^^')-i ^ 0(r - TS*i?^-^) holds even for rather small values of R (9 is the 



step function). Hence: 

PoAR)= / P{Y)dY. (6) 

Jo 

The overlap at the centre of the cavity is given by qdR) = Pm{R)(li + Pout{R)Qo, where qi 
is the self-overlap of the cavity (which is not necessarily 1 at finite temperatures). Using 
equation ([6]), we finally obtain 

POO 

q^{R)-q, = (q,-q,) dY P{Y). (7) 

At this point one needs to make some assumption about P{Y). An interesting parameteri- 
zation is: ^ 

W = ^(|)' exp[-{Y/Y^r] = -^exp[-{Y/Y^r] (8) 

{u > 0), which describes either a stiff distribution (small values of Y suppressed) at large 
z/ or a soft one (small values of Y enhanced) for small u. Y^ is the typical value of T. The 
standard mosaic picture corresponds to z/ — oo. Equation ([T]) now reads, explicitly: 

q,{R) -qo = (gi - qo) exp (9) 

where the usual RFOT relation between ^ and S is recovered: ^ = (Yc/TT,*Y^^'^^^\ As 
in the Adam-Gibbs treatement [28], the growth of the length is caused by the decreasing 
of the configurational entropy at low temperature. The expression ^ is the compressed 
exponential form ^ used above to fit our numerical data. Within this framework, the 
thermodynamic anomaly ( = v{d — 9) is directly related (at low temperatures) to the 
exponent v describing the surface tension distribution. 

A non trivial P{Y) could arise for two main reasons. One is strong pre-asymptotic 
corrections to standard RFOT: though the surface tension between two typical states could 
be self- averaging in the R oo limit, fluctuations can be state dependent and important 
at finite R. This is often what happens for interfaces in random media (see SI). The second 
possibility is that the effective interface tension fluctuates even in the large R limit. If the 
exponent 6 turns out to be less than — 1, as suggested by phenomenological arguments in 
ref. m then the interface should be a highly fluctuating object as e.g. in the droplet theory 
of spin- glasses [29]. We find that and therefore z/, increase when T decreases indicates that 
the crossover of qc{R) becomes sharper at larger sizes (see Fig. ([T]). This is compatible with 
a finite-size effect scenario. This behaviour is also expected within the RFOT scenario that 
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predicts a vanishing surface tension at the mode-couphng transition Tmc, which behaves as 
a spinodal point. Coherent amorphous order droplets should therefore be fractal around 
Tmc and compact below [30], which suggests an increase of the effective value of u as T 
decreases. A first principle RFOT computation of qdR) for the model we simulated would 
be very instrumental to clarify this issue. 



III. CONCLUSIONS AND OUTLOOK 



We have unveiled a qualitative difference between the high temperature and deeply super- 
cooled equilibrium regimes: the influence of boundary conditions propagates into the bulk 
on an increasingly large lengthscale upon cooling. Furthermore, the growth of this length 
is accompanied by a sharpening of the decay. We have developed a theoretical framework, 
based on a generalization of RFOT, that explains these results as a one-state to multi-state 
transition governed by the surface tension distribution. The sharpening of the decay at low 
temperature corresponds to more and more regions developing a large surface tension. From 
a more general perspective, our numerical results strongly support a mosaic picture where 
the super-cooled liquid is characterized by a "hidden" static order [32] on an increasingly 
larger scale upon cooling. Although this is a natural consequence of RFOT, other theoret- 
ical approaches may also account for these phenomena at least on a qualitative level: in 
particular the frustration limited domain theory fTS] and, perhaps surprisingly, some kinet- 
ically constrained models (see [18] for a discussion of this point). The main difference is the 
physical origin of the growing static lengthscale: within RFOT and at variance with other 
approaches, it is tightly linked to the decrease of the configurational entropy. Our work 
opens the way to a quantitative study of this issue and, hence, to a clear-cut test of RFOT 
as a valid theory of the glass transition. From a more theoretical point of view, several 
crucial questions remain elusive: can a RFOT-like transition exist outside mean-field? How 
precisely can amorphous metastable states be defined? Does the notion of effective surface 
tension between these states make sense? A definitive test of the mosaic scenario requires 
to find a way to measure directly this surface tension and its distribution in the deeply 
supercooled phase. Work in this direction is in progress. 
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IV. METHODS 
A. Simulation 

We have studied the soft-sphere binary mixture [23], a fragile model glass- former. In 
addition to capturing the essential features of fragile glasses, this model can be thermalized 
below the Mode Coupling temperature with the the swap Monte Carlo algorithm of [311 
Particles are of unit mass and belong to one of two species /z = 1,2, present in equal 
amounts and interacting via a potential 



where the radii are fixed by the conditions cr2/cri = 1.2, (2cri)'^ + 2(cri + cr2)^ + (2cr2)^ = 4£q, 
and io is the unit of length. The particle density is p = N/V = Iq^. A smooth long-range 



r > a, where a, Bij, and Cij are fixed by requiring continuity up to the second derivative of 
Vij{r). Temperature is measured in units of energy. To obtain the reference configurations, 
4 to 8 replicas of systems with = 2048 or N = 16384 where equilibrated in a cubic box 
with periodic boundary conditions. Then the overlap was computed in systems with frozen 
boundaries and M mobile particles, with M =20, 30, 40, 50, 100, 150, 200, 300, 400, 800, 
1600, 3200 mobile particles (1.684 < R < 9.142). The results are averaged over 8 to 32 outer 
states (reference configurations). Data were collected for at least 10 relaxation times (up to 
10^ Monte Carlo steps), after discarding an initial portion of at least one relaxations time. 
Error bars were obtained from a jacknife estimate from sample-to-sample fluctuations. 




cut-off is imposed setting Vij{r) = Bij{a — r)^ + Cij for r > Vc 



^/S and Vij{r) = Cij for 
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B. Overlap 



To obtain equation ([s]), write the partition function for the mobile cavity surrounded by 
pinning state a 

= e-f'^'f- + J2 e-^^'^--^^'^-, (10) 



SO that 

pi'^tiR) = ^ • (11) 

Introducing J\fa{f, T) = 6{f — f-y)S{T — T^-y), the sum can be written 

J2^-(3R''f,-PR^r^, = jdjIdTe-^^'f-^'^'^NU,^), (12) 

= JdfJdTe-^^'f-^^'™'''^%aiT\f), (13) 

where in the last equality we have defined Pa(T|/) = Afa{f,T)/Af{f), andA/'(/) = exp[S(/)] 
is the number of states with free energy /. Equation ^ follows approximating the integral 
over / with the saddle-point method, which picks / = /* as the dominant contribution. 
Pa{T\f) is the fraction of states 7 (inside the cavity) with free energy / and effective interface 
tension T, which we assume to be -R-independent [33] . 

Now the integral in equations ^ and Q can be simplified using the saddle point method, 
which is a very good approximation even for rather small values of R. Because of the 
exponential term in R^, the integral is dominated by the lowest values of T supported by the 
distribution pa{T\f). There are two possibles cases: (a) If this function has a left edge, i.e. it 
vanishes for T < ya, one finds up to subleading terms: / dT ]5q(T|/) e"^^^** ~ Q-PvaR" _ 
Otherwise: (b) there are arbitrarily small effective tensions T. Remarkably, in this case, 
depending on the form of Pa{T\f) at small T, one can obtain a renormalization of 6. For 
instance in the case Pq-(T|/) ^ exp(— Ca/T"), where a is a positive exponent, one finds 
/ dT Po(T|/) e'^'^^" ~ e^^^"-^" , where y'^ is a constant dependent of the outside state a 
and temperature and 9' = a9/{a + l) < 9. 

APPENDIX A: COMPARISON WITH THE RESULTS OF REF. 22 

Following the same numeric protocol of this work, in [22] the influence of the boundary 
conditions on the total overlap within the sphere was studied. It was observed that the decay 
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of '?tot(-R) is described sufficiently well in the single state framework and the standard RFOT 
scenario was ruled out. In this work we show however that at low enough temperatures 
neither the standard RFOT nor the single state scenario account for the local overlap data 
and a generalized RFOT theory is presented which is far more successful. 

One might retrospectively wonder if the data in [22] could have been used to discriminate 
between the one-state and the generalized RFOT scenarios. The answer is no. Allowing for 
a R dependence of go and qi of the type suggested in |22j: 



2(1 



x-^ 



+ %;1 



go;i(i?) = 3(l-g*;i) 
with X = R/Xq-i, the generalized RFOT prediction for gtot(-R) reads now: 
qtot{R) = qoiR) + (qiiR) - QoiR)) exp [- 



(Al) 



(A2) 



In the figure below we show that the generalized RFOT encoded in (A2) and the one- 
state prediction given by the formula (6) of ref. |22]fit the total overlap data at T = 0.203 
at a comparable level of accuracy. Following Occam's razor principle ( "entia non sunt 
multiplicanda praeter necessitatem") , one should choose one-state theory, which is the one 
with the smallest number of parameters. 
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